


AN EXPERIMENT TO MEASURE LASER BEAM WANDER 
AND BEAM SPREAD IN THE MARINE BOUNDARY LAYER 
NEAR SHORE 



David Al bert Bea 1 1 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 






AN EXPERIMENT TO MEASURE LASER BEAM WANDER 
AND BEAM SPREAD IN THE MARINE BOUNDARY LAYER 
NEAR SHORE 

by 

David Albert Beall 

Thesis Advisor: E. C. Crittenden, Jr. 

December 1973 



Approved public. xcJLcjl ie; dutnlbutlon unlimited. 



An Experiment to Measure Laser Beam Wander and Beam Spread 

In 

The Marine Boundary Layer Near Shore 



by 



David Albert Beall 

Lieutenant Commander, United States Navy 
'B.A., University of Washington, 1961 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN PHYSICS 



ru 




ABSTRACT 



A system to measure laser beam wander and beam spread in the atmosphere 
over the ocean has been designed, constructed, and tested. The apparatus 
employed a high resolution scanning telescope with a potential for use in 
a broad range of visual and infra-red wavelengths, and with the ability to 
measure beam wander and beam spread variations on the order of a few 
microradians in angle of incidence. Three successful trials with a pro- 
pagation path over. the southern end of Monterey Bay were conducted. Data 
was processed using analog and digital computers. RMS values of beam 
wander from 4.6 to 30.2 microradians were observed. 
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I. PROBLEM DESCRIPTION 



With the increasing knowledge and application of lasers in our society 
it has become important to understand more fully the ways in which a laser 
beam interacts with the medium through which it passes. Of particular 
importance are the effects of ordinary air, with its variations in density, 
temperature, humidity, and particulate content. There are two broad 
classes into which these effects may be categorized. The first class is 
only seen when high power beams are used, and are brought about when a 
significant amount of laser energy is transfered to the medium, causing 
meteorological phenomena over and above those that would be present normal- 
ly. These phenomena in turn have an effect on the beam. Thermal blooming 
and up-wind drift are examples of such non-linear disturbances. The 
second class are those effects which are due only to naturally occurring 
atmospheric fluctuations. Such effects may be observed in both low and 
high power beams. They include absorption, scattering. Modulation Transfer 
Function (MTF) variation, scintillation, beam wander, and beam spread. 

This thesis deals with the latter two phenomena as they occur in the 
atmosphere near the ocean surface. In particular, the development of an 
experiment to measure beam wander and beam spread in the atmosphere im- 
mediately above the sea surface will be described, together with some pre- 
liminary results of a small number of trials of such an experiment. 

It should be emphasized that the objectives of the work reported herein 
were somewhat limited from the outset. Due to the fact that the primary 
research ship, the R.V. ACANIA, was required for numerous other Naval 
Postgraduate School projects, that none of the apparatus was previously de- 
signed or assembled, and that the entire question of beam wander and beam 
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spread over the ocean was relatively unexplored, it was envisioned that the 
task of developing a fully comprehensive experiment, with the capability of 
operating over the open ocean and more than one optical wavelength, would 
require the efforts of two or more successive thesis students. Thus, this 
is just the first report on what is expected to be an ongoing project. 

A. MATHEMATICAL DESCRIPTION OF BEAM WANDER AND BEAM SPREAD 

Since the index of refraction of a medium depends upon its density, a 
beam passing through a medium with variations in density may be expected to 
follow a path which is not straight. If the diameter of the beam is not 
much less than the transverse distance in which significant variations 
occur, then the bending of the beam will be fairly complicated, since 
different portions of the beam will pass through regions of different re- 
fractive index. Thus, the beam will change shape, change direction, and 
exhibit regions within its cross-section where interference causes greater 
or less intensity than would be expected in a medium with constant density. 
These effects are known respectively as beam spread, beam wander and scintil- 
lation. Obviously, these same effects will also occur in any beam regard- 
less of diameter if significant variation in density along the path of 
propagation is present. In the atmosphere normal meteorological conditions 
commonly cause both longitudinal and transverse density variations over a 
range of scales sufficient to induce such effects in beams of all sizes. 

In fact turbulence scales near the earth's surface may be expected to range 
from the order of one millimeter to a few meters [Ref. 1]. The term, 
turbulence scale, is somewhat loosely defined. It may be roughly visualized 
as the diameter of a turbule, or region where the air rotates about a 

3 

common point. Over a path of, say, 1km, then one might expect at least 10 
to 10^ such turbules. 
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Clearly a statistical model is the only formulation which could be 
expected to describe the effects of turbulence on a laser beam. The only 
successful models existent require the assumption of randomness of trans- 
verse refractive index fluctuation. The preeminent basis for such models 
was introduced by Tatarski [Ref. 2]. First, C n , the refractive index 
structure constant, is defined. 

C n (f) = [0 n (?)] 1/2 r- ,/3 , l 0 <r<L 0 , (1) 

where 

_ D n (r) - <[n(R)-n(R)+F)] 2 > av . (2) 

Here, r = a vector between two points in the transmitter aperature, 

R = a vector to any point within the conmon area of two circles 
with centers at the endpoints of r, 

n = the index of refraction at a point, 

1q = the smallest turbulence scale present, 

l_ 0 = the largest turbulence scale present, 

av = indicates an ensemble average over all realizations of the 
atmosphere for fixed meteorological conditions. 

Fortunately C n can be measured, so that analytical solution of (1) and 

(2) is not required. This can be done either through measurement of 

scintillation or micrometeorological observations. Recent measurements of 

C n by Haagensen [Ref. 3] demonstrate that meteorological and optical 

measurement of C over the ocean surface are well correlated, 
n 

Tatarski showed that a power spectrum based on refractive index fluc- 
tuations is proportional to the velocity spectrum of the turbulence [Ref. 2], 
From this fact, plus the assumption of isotropy which allows the conversion 
of (1) and (2) to scalar equations, Strohbehn developed an expression for 
a refractive index power spectrum, 4> n (K), [Ref. 1]: 
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(3) 



♦ n 00 = 0.033C n ^exp(-K 2 /K rn 2 )K 0 ~^ /3 (l+K 2 / K Q 2 )~ 11/6 

where, 

K = spatial wave number, 

•S„ = 5 - 92 'V 

K q = 2-/L 0 . 

From 4> n (K), a wave structure function, D(r,z), can be developed 
[Refs. 1, 2, 4], 

D(r,z) = 4irk 0 2 z/Jdy/dK<{i n (K){l-exp[iKr(l-y)]}, 

where, 

z = distance along the beam from the transmitter, 
kg = optical frequency wave number. 

D(r,z) can be inserted in the form of a weight function into the 

Fresnel -Kirchhoff integral to give a time averaged irradiance in the re- 

ceiver plane, 

I(w) = (kQ 2 /2Trz)/drF(r,z)exp[-l/2D(r,z)] exp[-iKw*r/z] (5) 

where, 

w = position in the receiver plane, 

F(r,z) = the aperature function projected through the distance z. 

That is, 

F(r,z) = /dRU(R+l/2r)U*(R-l/2r)exp(ikQR-r/z) (6) 

where U is the amplitude of the radiation in the aperature [Ref. 4]. 

As pointed out in Ref. 4, equation (5) is not susceptible in general to 
analytic solution for beam wander and spread. However, several approxi- 
mations have been made [Refs. 4, 5, 6, 7, 8, 9]. These generally give an 

expression for beam wander or beam spread in terms of C n> z, and either Kq 

or 1 q. Dowling and Livingston have tabulated many of these expressions in 
common units and dimensions in Ref. 4. For example, their own expression, 
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based partly on work by Fried, gives beam wander, defined as the mean 
square motion of the energy centroid as 

* 2 = 32ir 2 zf/3L 1 , (7) 

where, 

L ] = 1/[1/2X0. 0330C n 2 (l Q /5. 92) 1/2 ], (8) 

f = a complicated function of transmitter aperature and 

turbulence developed by Fried, which is independent 
of z, k Q , and C n [Refs. 4, 10]. 

Dowling and Livingston also show that assuming a gaussian beam, 

e' 2 = e 2 +«j> 2 , (9) 

where, 

6 = short term average beam width expressed as the angle 
in radians from the transmitter subtended by the 1/e^ 
power distribution width in the focal plane (waist) 
of the beam, 

Q' = long term average beam width expressed in tftie same 
terms [Ref. 4]. 

In figure 1 of Ref. 8, Chiba gives a partly qualitative graph, which 
would seem general enough to cover all the approximations of the Tatarski 
model in representing the increase of beam wander witto turbulence. 

B. BROAD DESIGN CRITERIA FOR THE APPARATUS 
1 . Uniqueness of the Experiment 

At the outset of the investigation no other experiment to measure 
beam wander or spread over the open ocean was known to have been conducted 
or to be in progress. This placed certain constraints on design. In 
particular, the relatively high probability that any particular complete 
experiment would not be successful indicated a fairly cautious approach. 
Each component would be built, tested indoors, and possibly rebuilt and 
retested before inclusion in the experimental system. The entire system 
was tested in a 400 foot hallway before the first over-water trial. 
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2. Portability 

A fairly compact portable apparatus was required for three reasons. 
First, as mentioned above, it was deemed important to test the apparatus 
indoors preliminary to over-water experiments. Second, although to date 
the experiment has been conducted successfully across bays, that is with 
the apparatus on one shore and the transmitting laser on the other, it is 
desirable that it ultimately be conducted between two ships if possible 
or at least from ship to shore. This is due to uncertainties as to the 
effect some shorel ine 'profiles may have on atmosphere flow patterns. Third, 
it was not known whether the relative unchangeability of weather conditions 
in the vicinity of Monterey would provide a sufficient base for correlations 
with beam wander and spread. Thus, it might have been necessary to conduct 
the experiment in other locales. 

3. Precision 

Overland measurements of beam wander by Dowling and Livingston 
[Ref. 4], and by Chiba [Ref. 8] indicated that variations in beam direction 
on the order of 100 microradians for path lengths in the kilometer range 
might be expected. This indicated the requirement for high quality optics, 
precision scanning, and a means for compensating for motion involved in 
tracking a moving target, i.e., the shipboard laser. In fact, RMS beam 
wander measurements on the order of 10 microradians have been generally 
observed on cross bay experiments. 

4. Scan Repetition Rate 

Concerning the possible characteristic frequencies which might be 
associated with beam wander, data from Ref. 4, as well as subjective 
observations of laser beam behavior on previous experiments, such as the 
scintillation work by Haagensen [Ref. 3], indicated that frequencies in the 
1 to 10 Hz regime might be expected. Thus at a minimum, the apparatus 
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should be capable of scanning the beam or an image of the beam on the 
order of 100 times a second. 

5. Wavelength Capability 

The experiment as designed and conducted to date has utilized a 
source of 6328 Angstrom He-Ne laser radiation. However, it is desirable 
that in the future the apparatus be used in the infra-red range out to 10.6 
microns. Thus it was designed as much as possible with front surface re- 
flective optics to provide a broad band capability. 
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II. DEVELOPMENT OF THE APPARATUS 



A. BASIC CONCEPT 

As noted above, beam wander and spread measurements over land were made 
by Dowling and Livingston in 1972 [Ref. 4]. Their apparatus was straight- 
forward, consisting of two laser beams, one at 6328 Angstroms and one at 
10.6 microns, transmitted by a Newtonian telescope to a 91.5 cm diameter 
spherical mirror at the receiver site. The receiver was located at the 
focal plane of the beam. The radiation reflected from the 91.5 cm mirror 
was brought to a small diameter on a galvanometer mirror, which in turn 
sent images to two slits, one for each wavelength, behind which detectors 
were located. In effect an image of a small portion of the 91.5 cm mirror 
was placed at each slit. With the galvanometer mirror movement driven by a 
triangle waveform current, an entire diameter of the mirror was scanned at 
a repetition rate of 75 Hz. One of the two systems reported by Chiba 
[Ref. 8] employed a position sensitive detector. The other was largely 
similar to the system used by Dowling and Livingston. 

The relative lack of portability, questionable ability for tracking a 
moving target, and cost of a high quality concave mirror of the size re- 
quired obviated the above approach for an early over-water experiment. 

Rather it was decided to use the principle of reversibility in the Naval 
Postgraduate School apparatus. The telescope would be placed at the re- 
ceiver, rather than the transmitter, and the image to be scanned would be 
the image of the source, rather than that of the beam focal plane. Assuming 
that the turbulence is symmetric in the propagation direction about the 
mid-point of the propagation path, the principle of reversibility yields 

that the motion of the image of a stationary laser cavity due to atmospheric 
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refractions is the same as if the image were stationary and the object in 
motion. The diffraction limited image gathered from defocused light of the 
laser is equivalent to the intensity distribution of the focused laser 
radiation when viewed at its waist, i.e., its diffraction limited focal 
plane. 

Thus the basic concept was to use a divergent laser beam as a source, 
and record the motion and shape of the image of the laser crossover point 
(Figure 1) formed by the telescope. 

B. EARLY CONFIGURATION OF APPARATUS 

1 . Telescope 

The first telescope built used a 9.2 cm diameter, 92.4 cm focal 
length parabolic mirror previously used by Hildebrand [Ref. 11] in atmos- 
pheric MTF research. This was combined with a 3.8 cm diameter, 24.5 cm 
focal length convex spherical mirror and a flat 90 degree turning mirror in 
a Cassegranian telescope. The arrangement is shown in Figure 2. This 
telescope had an equivalent focal length of approximately 6.6 meters. 

2. Scanner 

The first scanning device built employed an eight sided prism which 
was rotated at about 60 Hz. A magnetic field was provided perpendicular to 
the axis of rotation. The prism mounting had eight narrow grooves parallel 
to the rotation axis, one beneath each corner of the prism. A tape recorder 
magnetic playback head was located adjacent to the prism mounting. As the 
grooves on the mounting passed the playback head, they caused perturbations 
in the magnetic field, which activated a signal in the playback head. This 
signal was used to trigger an oscilliscope. The prism was positioned in the 
output optical path of the telescope. The rotation of the prism caused the 
image formed by the telescope and prism to be scanned across a slit, behind 
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FIGURE 1. Transmitter Optics and Propagation Path 
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FIGURE 3. Eight Sided Prism Scanner 
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which a detector was mounted. The signal from the detector was fed to the 
vertical input of the oscilloscope and to a tape recorder for storage. 

Thus, if the grooves in the prism mounting could be located precisely enough, 
and the rotation speed could be kept constant, consecutive images scanned 
past the slit would appear as optical intensity curves on the oscilloscope 
referred to a fixed position in the image plane and be recorded. Figure 3 
shows this arrangement. 

3. Advantages and Difficulties 

The practical problem of cutting the grooves in the prism mount and 
positioning the prism with sufficient accuracy to stabilize the image on 
the oscilloscope proved insurmountable. Therefore it was necessary to use 
only one out of every eight images transmitted, in order that the image 
displayed have a constant phase relationship with the magnetic groove trigger 
signal. Thus only about eight usable images per second were available, 
obviously too few for good resolution of a phenomena expected to have char- 
acteristic frequencies in the 1 to 10 Hz range. In addition, the prism 
obviated the capability to operate in the visual and infra-red regions simul- 
taneously. No material with good transmission and suitable mechanical pro- 
perties in such a broad band was known to be available. 

The secondary mirror was manufactured by silvering a lens to provide 
near 100% reflectance. This mirror produced some astigmatism in the image 
in certain directions. Although good quality convex mirrors are available, 
the other considerations mentioned above caused a reconfiguration of the 
apparatus, and the Cassegranian telescope was discarded in the process. 

Despite these faults, this configuration had many strong features. 

As foreseen by Hildebrand in Ref. 4, it eliminated the nuisance of having 
successive images displayed with reversed positive direction of the abscissa, 
which occurs in galvanometer type scanners. The Cassegranian telescope 
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allowed a relatively high magnification in a reasonably short telescope, 
i.e., the ratio of the focal length to the telescope length was about four. 
Despite the problem of astigmatism, the magnification achieved produced 
excellent structure resolution in the image. Figure 4 shows some of the 
images produced by this system. The object was the optical crossover point 
of a 6328 Angstrom, 0.1 mW TM 00 laser with divergent optics which provided 
a virtual object distance of 104 m. The laser was located about four fifths 
of this distance from the telescope in a hallway with numerous ventilators 
in the overhead. Scintillation and beam wander are clearly visible. Ex- 
posure time was greater than the scan time, so that more than one image is 
seen in each photograph. The less bright lines are persistent traces of 
images displayed before the camera shutter was opened. 

C. PRESENT CONFIGURATION OF APPARATUS 

Despite the advantages of the system described above, its disadvantages 
were too great to be overcome. Thus, a new system has been developed, which 
employs galvanometer mirror scanning and a Newtonian telescope. The 
primary mirror is a 25.4 cm diameter 204.5 cm focal length spherical mirror. 
Its reflected light is intercepted just short of the focal plane by a flat 
mirror oriented 45 degrees from the telescope axis, and transmitted through 
a filter and diverging lens at the top of the telescope. Here it is turned 
by another 45 degree mirror and brought to a focus in front of the scanning 
mirror. It is reflected from the concave scanning mirror, passes back be- 
tween the upper folding mirror and the diverging lens, and is brought to a 
focus on a vertical slit. A detector is mounted directly behind the slit. 
The positions of both folding mirrors, the diverging lens, the galvanometer 
mirror, and the slit/detector assembly are adjustable so that an image may 
be formed at the plane of the slit for object distances from 100 meters to 



20 





Fir.i'PF 



’'sares Obtained 1 r a 

And f i c‘ p - J r-r i 

b'ori ront.o ’ I'rs ' j c. ■ r 



seer? rian T e’ t-scoiy 
>■ T r. * a r s i * '■ vs. 
c ! con Hary-. 



dl 



Infinity. In addition, the mounting for the upper folding mirror, scanner, 
and slit/detector assembly can be rotated through a complete circle in the 
horizontal plane. This allows the beam to be scanned through a linear path 
at any angle in its transverse plane without rotating the whole telescope. 
Figures 5, 6, and 7 show the telescope and scanner assembly. 

To operate the galvanometer as a scanner, it was necessary to replace 
its suspension with a stiffer wire, since the original gold foil suspension 
was not expected to be rugged enough to survive transportation to the field. 

The sinusoidal drive current causes the galvanometer movement to act as 
a forced harmonic oscillator, with the input current as the forcing function. 
Thus the angular deflection of the mirror is sinusoidal in time, with the 
same frequency as the drive current. The phase angle between the driving 
current and the mirror deflection is constant once steady state has been 
achieved. [Ref . 12]. Tests of the scan mirror showed that the center of the 
angular deflection remained constant to within 0.004 of the amplitude of 
the deflection. For the instrument settings used in subsequent over-water 
experiments, this corresponds to an accuracy of about 2.8 microradians in 
beam wander. 

• To provide fine steering of the galvanometer, a variable D.C. bias 
voltage is applied to the driving current. This has the effect of shifting 
the center of the scan. It was envisioned that this bias voltage could be 
removed during data reduction. However, up to the present, this voltage has 
been used only to set the scan center prior to starting data collection, and 
variation has not been required during collection itself. 

The detector is an E.G.&G. Incorporated SGD-444 silicon photodiode. The 
slit opening is adjustable and has been used at various settings during the 
experiments conducted. The largest slit opening used has been 0.066 mm, and 
the preponderance of the data has been collected with a slit opening of 
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0.045 mm. Slit calibration was obtained with a traveling microscope. The 
length of the slit is 1.5 cm. 

Since there is no provision for eliminating motion of the image along 
the length of the- slit due to beam wander and spread in the direction per- 
pendicular to that being scanned, the question of the uniformity of the 
response of the detector is important. This has been measured with the 
result that a constant intensity beam trained along the length of the slit 
produced a maximum deviation of 2.5 percent in the detector response. This 
figure may be somewhat low, however, since the beam used was about twice the 
diameter of the images commonly seen in the field experiments, thus allowing 
some averaging effects. 

The entire scanning/detecting assembly is provided with a cover to re- 
duce side and back light. The cover is not completely light tight, however 
this is only a matter of convenience, since low frequency intensity levels 
are removed during data processing. 

The detector signal is displayed on an oscilloscope, which also amplifies 
it and outputs it to a bias control circuit. It is then passed through a 
10:1 amplifier to a Hewlett-Packard 3960E F.M. four channel instrument tape 
recorder. The recorder is operated at 15 inches per second, providing 0 to 
5000 Hz frequency response. Thus, it is possible to record 10,000 intensity 
levels per second. The recorder responds to input voltages in the range 
-10 to 10 volts. Since the image intensity is always a positive signal, the 
bias control circuit allows use of the full voltage range of the recorder. 
This is done by adjusting the bias so that the D.C. level of the input 
signal is at the minimum (-100% on the recorder level indicator) of the 
voltage response range. This procedure effectively increases the signal to 
noise ratio of the recorded signal. 
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The scan drive signal is also recorded, and at the same time used to 
trigger the oscilloscope. Figure 8 shows a schematic diagram of the entire 
system. 

Calibration of the system is accomplished using a coarse diffraction 
grating. To obtain a calibration for a particular set of adjustments on 
the telescope, the grid is placed in the incoming beam directly in front of 
the telescope. The grid bar center-to-center separation is 8.95 mm. Thus 
the distance between the diffraction peaks on the oscilloscope, when using 
6328 Angstrom light, corresponds to 70.7 microradians angular deflection of 
rays entering the telescope. If the X axis of the oscilloscope is driven 
by the scan drive signal, then the diffraction pattern appears on a linear 
horizontal scale, and may be used as a calibration. Figure 9 shows the 
diffraction grid and a typical diffraction pattern it produced. 

All of the trials of the experiment discussed below employed the config- 
uration of apparatus just described. 
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III. EXPERIMENTS 



To date, five over-water experiments have been attempted. Only three 
of these produced successful results. The other two, San Micholas Island, 
September 1973 and Monterey Bay, November 1973, are discussed briefly, be- 
cause it is felt that the difficulties encountered can be avoided if known 
beforehand. 

A. SAN NICHOLAS ISLAND, SEPTEMBER 1973 

This experiment was originally scheduled as only a scintillation experi- 
ment, with the objectives of testing various details of scintillation 
theory [Ref. 13] and initially testing a gyro stabilized ship mounted laser. 
The beam wander and spread experiment was added when it became apparent 
that the apparatus would be ready in time. Two transmitters were available, 
one the ship mounted laser aboard the R. V. ACANIA, and the other located 
ashore across Laser Bay from the receiver site, located on the northwest tip 
of San Nicholas Island in the Santa Barbara Channel. Equipment was trans- 
ported to the island by air, and was set up at the site in tents. 

Despite all efforts to protect the equipment, a northwest wind of up to 
30 knots continually blew sand and salt spray into the tents, making it 
necessary to keep the apparatus wrapped in plastic and non-functional most 
of the time. At times, visibility was reduced to a few hundred meters by 
the amount of sand in the air. In addition, the size of the tents required 
the front end of the telescope to protrude from the tent door. During 
attempts to take data the wind caused so much vibration of the telescope that 
stabilizing the scan was impossible. These conditions, combined with crew 
fatigue and unforeseen logistic problems, rendered the situation untenable, 
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and the experiment was cancelled. Lessons learned from this trial were as 
follows. 

The apparatus should be enclosed in a solid structure. For future 
trips to San Nicholas Island a van should be used. 

The shore based crew should arrive and set up well in advance of the 
ship's arrival. At least two days advance time to set up and test equip- 
ment and become acclimatized to the harsh weather conditions on the island 
should be planned. 

Any feasible measure to provide personnel support should be taken. The 
equipment van should be provided with cooking and sleeping facilities, 
since the nearest messing and berthing facilities are a 45 minute drive 
from the receiver site. The only transportation is an open army weapons 
carrier, and much of the trip to the air terminal/messing/berthing area is 
over dirt roads that are in extremely poor repair. This fact should be 
kept in mind when preparing equipment for transportation. 

There is no pier on the northern part of the island, and beaching a 
whale boat is hazardous, since the bay is only slightly protected. Two 
personnel exchanges between the ship and shore were accomplished, but they 
required wading out chest deep to the ship's boat. A zodiac boat would 
probably be capable of beaching. 

B. MONTEREY BAY, NOVEMBER 1973 

Throughout the development of the beam wander and spread receiver ap- 
paratus, a parallel project to manufacture a gyro stabilized laser trans- 
mitter for shipboard use was being carried out under Professor Eugene C. 
Crittenden, Jr. One of the intermediate tests of the gyro system was con- 
ducted in November 1973 with the beam wander and spread apparatus as one of 
the receivers used. A manual tracking system was attached to the telescope. 

The stabilized transmitter was located aboard the R. V. ACANIA at a distance 
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from the receiver site of about 1.0 km. Although the stabilization system 
operated properly, vertical motion of the source due to swells, two to four 
feet, as well as lateral ship motions were too great to keep the image in 
the field of view of the telescope. 

Based on this attempt, it appears that fully automatic tracking by the 
receiver will be required for future ship to shore experiments. 

C. MONTEREY BAY, 9, 18, AND 22 OCTOBER 1973 

These three successful trials were conducted with the transmitter located 
at the Holiday Inn Motel, in the northern outskirts of Monterey, and the 
receiver at Stanford University's Hopkins Marine Station, Monterey. The 
range is shown in Figure 10. The source was a 6328 Angstrom TM 00 laser 
with 3 mW output. The propagation path was horizontal and was 4.1 km in 
length. Transmitter and receiver were located 15 feet above the mean high 
water mark. 

On the 9th and 22nd, no meteorological data over the propagation path 
was available. On the 18th, meteorological C n measurements were taken 
aboard the R.. V. ACANIA stationed near the mid-point of the propagation path. 
However, processing of that data has not been completed. Optical C n measure- 
ments were taken for all three days, and have been processed [Ref. 13]. 

On the 9th and 18th, the prevailing wind was from the west at about 10 
knots. Cloud cover varied from 0 to 100 percent. On the 22nd, data was 
taken after a morning shower. The atmosphere was exceptionally clear. The 
prevailing wind was from the south at approximately 20 knots, and the sky 
was clear throughout the trial. 

Data was taken quasi-continuously on the 9th and 18th. On the 22nd, 
data was taken in 50 second periods at 10 minute intervals. The experiments 
were conducted from 1754 to 2314 9 October, from 1855 to 2154 18 October, 
and from 1555 to 1754 22 October. 
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FIGURE 10. Chart of Monterey Ray Test Range 
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A scan frequency of 120 images per second was used on all three days. 
As mentioned above, slit openings varied from 0.045 mm to 0.066 mm, with 
most of the data taken with the former value. 
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IV. DATA REDUCTION 



The data is processed in five major steps, which are described below. 

A. ANALOG TO DIGITAL CONVERSION 

The tape recorded optical intensity and scan drive signals are converted 
to digital information in the Naval Postgraduate School Electrical Engineer- 
ing Department Computer Laboratory, using a COMCOR 5000 analog computer and 
a Scientific Data Systems SDS 9300 computer. When digitized the data are 
stored as four byte integers in octal representation on seven track magnetic 
tapes. Operating procedures for the converter and computer may be found in 
Refs. 14 and 15. 

Several programs for A/D conversion are available, and any facility which 
uses the SDS 9300 will no doubt have such programs available. They are 
written in XDS FORTRAN, a variant FORTRAN language, and make extensive use 
of assembler language subroutines. For these reasons, they are not included 
in this thesis. However, for convenience of Naval Postgraduate School users 
definitions of variables which must be input by the operator are given below. 

NCHAN is the number of channels to be recorded. For this experiment 
NCHAN is two, however, the A/D system is capable of handling up to eight 
channels. 

NREC is an upper bound on the number of records that will be written on 
the tape. A record is a block of data recorded on the tape separated from 
the previous and subsequent record by a space. Division of the data stream 
into records allows reading and writing on tape to proceed quickly. 

NSAMP is the number of samples of each channel per record. NSAMP equal 
512 has been used with good results. 
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ITAPE is the magnetic tape unit number which is set on the magnetic tape 
console when mounting the tape. 

NF is the number of files on the tape to be skipped before commencing 
recording. A file is a complete set of converted data, i.e., all the data 
from the interval over which beam wander and spread are to be computed. An 
end of file (EOF) mark is written on the tape at the conclusion of each 
file. The EOF mark cannot be passed by a computer when reading the tape, 
and will trigger the END option of a FORTRAN read statement. 

The maximum sample rate at which the system can operate is dependent on 
the number of channels being digitized. For two channels a sample rate of 
3333 samples per second on each channel is possible. However, better re- 
sults in subsequent processing have been obtained if a sample rate of 2500 
is used. Digitizing is done with the original recording played back at 1/4 
speed, i.e., 3 3/4 inches per second. Thus, each channel is sampled at 
10,000 per second real-time, that is the time reference frame in which the 
data is collected. This corresponds to the frequency response of the tape 
recorder. This is believed to be the reason why better results are obtained 
at the lower sampling rate. This sample rate allows about five files of 
15 seconds real-time to be recorded on the 2400 foot seven track computer 
tape. 

B. OCTAL TO HEXADECIMAL CONVERSION 

It is possible to conduct all subsequent data processing steps with the 
SDS 9300 computer. However, it is understandably the policy of the computer 
laboratory only to undertake jobs, other than experimental computer science 
work, which cannot be accomplished by the school's general purpose computer, 
IBM 360/67 system. Thus, it is necessary to convert the octal integers 
output by the SDS 9300 to hexidecimal representation for use by the IBM 360. 
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This is done using the FORTRAN program 7T09 listed below. The program is 
executed in the IBM 360 computer. It is based on a similar program given 
in Ref. 16. 

Logical unit numbers used in all FORTRAN READ and WRITE statements in 
this thesis are as follows: 

2: seven track tape 

4: nine track tape 

5: card reader 

6: Line printer 

7: card punch. 

The program converts the octal integers to hexadecimal integers, con- 
verts these to four byte real numbers, and stores them on a nine track 
magnetic tape. The actual numerical conversion is accomplished by the 
assembler language subroutine, FORM, which has been precompiled and resides 
in the user library, MPSLIB [Ref. 16]. 

It would be possible to conduct all subsequent data processing steps 
without re-writing the data on the nine track tape. However, this would 
approximately double the computer time in operating with the seven track 
tape. Since there is only one seven track mounting in the IBM 360 system, 
and several other users, it is advisable to go through the intermediate step 
of writing the nine track tape. This also facilitates an earlier receipt of 
portions of the final data processing product, since the nine track tape 
data can be conveniently processed one file at a time. 

The program 7T09 requires approximately seven minutes CPU time to pro- 
cess five files with NCHAN equal two and NSAMP equal 512, where the files 
contain 15 seconds of real-time data. 

Significant variables in the program are defined below. 
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NRECL is the total number of samples per record. It must be equal to 
NSAMP times NCHAN. 

FACTOR is an arbitrary constant introduced to provide convenient sized 

values of the intensity and drive signals for later processing. With the 

instrument tape recorder output in the range of 1-10 volts during A/D con- 

-4 

version, FACTOR equal to 10.0 produces data values in the range 10.0 to 
1000.0, depending on the amplifier settings that are used on the C0MC0R 5000 
computer. FACTOR equal to 100.0X2.0 will yield the actual voltage at 
the A/D trunk of the COMCOR 5000. 

The arrays IDAT and DAT must be dimensioned to NRECL. 

The FORMAT statement 15 must also provide for reading exactly NRECL 
words, i.e., as written in Appendix I the statement calls for reading 

8 X 128 = 1024 = NRECL 

samples. 

C. COMPUTATION OF BEAM WANDER AND IMAGE SHAPE 

Listed below is the FORTRAN program BEAM WANDER/SPREAD, which computes 
RMS beam wander, time average image shape, and time average image shape with 
beam wander removed. The main sequences of the program are given in 
Figure 11 and are shown in detail on Figure 12. 

The program consists basically of two nested loops. The outer loop 
performs operations on each record, and the inner loop operates on each 
image within the record. The beginning of the inner loop is signified by 
the box marked branch in Figure 11. 

After initializing arrays and variables, the program conmences the outer 
loop by reading the first record from the magnetic tape. It sorts the data 
so that the image intensity signal is placed in array Cl (for channel 1), 
and the scan drive signal in channel 2, C2. It then searches for the first 
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Initialization 







FIGURE 11. Beam Wander/Spread Master Flow Diagram 
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Initialization and Starting Main Loop 




FIGURE 12. Detailed Flow Diagram 
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Reading and Sorting Data 




FIGURE 12. Continued 
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Finding First Maximum of Channel 2 




FIGURE 12. Continued 
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Finding First Minimum of Channel 2 
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FIGURE 12. Continued 
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Removing Bias From Channel 2 
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FIGURE 12. Continued 



Finding First Zero of Channel 2 




FIGURE 12. Continued 
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FIGURE 12. Continued 
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Normal Case: Removing Bias From Channel 1 
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FIGURE 12. Continued 
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Normal Case: Removing Bias From Channel 1 (Cont.) 




FIGURE 12. Continued 
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Normal Case: Finding Period of Channel 2 







FIGURE 12. 
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Continued 



Normal Case Integrations 





FIGURE 12. Continued 
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Overflow Case: Loading Overflow 
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Overflow Case: Removing Bias From Channel 1 
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FIGURE 12. Continued 
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Overflow Case: Finding Period of Channel 2 
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Overflow Case: Finding Period of Channel 2 (Cont.) 




FIGURE 12. Continued 
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Overflow Case: Integrations 




FIGURE 12. Continued 
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Overflow Case: Loading Arrays and Preparing Next 

Integration Sequence 
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Loading Overflow at End of Record 
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FIGURE 12. Continued 
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Output 
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maximum and minimum of C2, and if it finds them it computes the bias of C2 
and subtracts it off, leaving a purely A.C. sinusoid. If the program cannot 
find the maximum or minimum of C2 due to abnormalities in the sine wave, 
it returns to the start of the loop and reads the next record. Normally 
this condition doesn't occur, however, and the program continues in the 
outer loop by finding the first zero of C2. It does this by searching for 
the smallest absolute value of C2 in a domain slightly larger than a half 
period of the sinusoid as determined from the points where the first maxi- 
mum and minimum appeared. Again, this is normally successful; however, 
if it is not, the program advances another half period and searches again. 

At the branch point, the program determines whether there are any samples 
in the overflow arrays OFLOl and 0FL02. These arrays are filled at the end 
of a record if the domain of Cl required for the integrations described 
below overlaps the end of the record. At the beginning of a record, if 
there are values in these arrays, the program completes the loading of them 
and conducts the inner loop operations on the first image using 0LF01 and 
0FL02 instead of Cl and C2. All subsequent images in the record are pro- 
cessed using the values in Cl and C2. The inner loop operations are as 
follows. First the bias of Cl is removed over a domain immediately before 
the domain of integration. Both of these domains have length specified by 
the user during initialization. The center of the domain of integration is 
always the same location as the zero of C2 found previously. After removing 
the bias from Cl over the domain of integration, the program searches for 
the next zero of C2. The difference between the current zero and the next 
zero is taken as a more accurate value of the half period, since C2 is a 
more rapidly changing function near its zeros than near its maximum or 
minimum. This value of the half period is then used in approximating the 
derivative of C2 for use in the integration sequence. The integration 
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sequence, which computes the centroid of the image intensity relative to 
the zero of C2, is described in detail in section 2 below. After this 
computation the program stores the position of the centroid in the array 
WANDER. The position of the centroid is then subtracted from the position 
of each intensity sample in the image and the samples are added to the 
array SPREAD. Each image is added to SPREAD with its centroid at the mid- 
point of the array. Thus, SPREAD accumulates a time average image shape 
with beam wander removed. In a similar way, but without subtracting the 
position of the controid, the program adds the intensity samples from the 
image to the array UNFXD to accumulate a time average image shape with beam 
wander included. Having completed these operations, the program advances 
to the next zero of C2 and returns to the beginning of the inner loop. It 
checks whether there is sufficient room left on the record to recycle through 
the loop. If there isn't, it checks whether the end of the record falls 
within the integration region of the next image, as described above. If not, 
the program returns to the beginning of the outer loop immediately. How- 
ever, if the end of the record falls within the integration domain, the 
program loads the overflow arrays from the beginning of the bias removal 
domain to the end of the record and then returns to the beginning of the 
outer loop. 

When the end of the file is reached, the program branches to the output 
sequence. Here it outputs the values stored in WANDER, SPREAD, and UNFXD 
to the line printer as well as to punched cards for later processing. 

If desired, another loop can be added encircling both the loops de- 
scribed above to allow processing of more than one file at a time. This has 
been done successfully with some of the data processed, but is not included 
here in order not to unnecessarily complicate the description. The start 
and stop points for such a loop are immediately after the first REWIND 
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statement and before the last REWIND statement, respectively. If used, this 
procedure requires DO loops to re-initial ize SPREAD, WANDER, and UNFXD at 
the beginning of each file. 

1. Initial ization 

The initialization sequence dimensions arrays, reads input parameters, 
rewinds the tape to the beginning of the file specified by the Job Control 
Language (JCL) statements submitted with the job [Ref. 17], and initializes 
the counting parameters defined below. 

L is the current record number being processed. 

M is the sequential number of the image being processed. 

IOFLO is the number of data points loaded into overflow arrays 0FL01 
and 0FL02 at the end of a record to provide continuity at the beginning of 
the next record. 

IERNUM is the number of times the program has failed to find a zero, 
maximum or minimum of the scan drive signal (channel 2). 

Variables which must be input via the card reader are as follows. 

IAP2PI is the approximate number of samples in one period of the 
scan drive signal. This may be estimated either from the sample rate used 
when digitizing the data or from the printed output of 7T09. 

NCHAN and NSAMP are as defined above. 

MAXERR is an upper bound on IERNUM. If IERNUM reaches this value 
the program automatically goes to the output sequence. 

INTREO is the half width of the integration domain to be used in 
computing the centroid of each image, expressed as a number of samples. Its 
value can be chosen by refering either to the 7T09 printout or a dual trace 
display of both channels to determine an upper bound on the fraction of the 
period of channel 2 required for the image signal to rise from and return to 
its D.C. level. The integration is carried out from Ig-INTREQ to Iq+INTREQ, 
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where Iq is such that C2(Ig) is nearest to zero for each image. In this 
regard it should be noted that the scan drive signal is phase shifted during 
the digitizing process to correspond to the phase of the deflection of the 
scanning mirror, so that images are entered in the data stream near each 
zero of channel 2. 

IBIREQ is the number of samples of the intensity (channel 1) required 
to determine the bias of channel 1 at each image. 

IAGREQ is the allowable variation in the length of half a period of 
channel 2, expressed as a number of samples. 

KO is the mid-point of the arrays used to accumulate time average 
image shape and time average image shape with beam wander removed. 

All arrays provided are one dimensional. Their functions and 
dimensions required are described below. 

DAT stores all the data from each record. It must have exactly the 
same dimension as was used for DAT in 7T09. 

SPREAD and UNFXD store the time average image shape with wander 
removed and time average image shape without wander removed respectively 
as they accumulate. The points in SPREAD and UNFXD are in the spatial 
rather than the time reference frame. That is, they are the nearest integer 
to the value of channel 2. Thus, they should be dimensioned to correspond 
to whatever fraction of the amplitude of channel 2 will contain the entir- 
ity of each image. For example, suppose either from observing the original 
tape recording or the printed output of 7T09, it is determined that every 
image occurs within tt/ 8 of the zeroes of channel 2, and that from the 7T09 
printout it is determined that channel 2 varies from -100 to 100. Then 
SPREAD and UNFXD should be dimensioned to at least 2 x 100simr/8. In general 
it is possible to dimension these arrays larger than one would ever need, 
and not worry about them thereafter, since they are small arrays anyv/ay. 
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WANDER stores the position of the centroid of each image, and should 
be dimensioned as large as the number of images in the file. 

Cl and C2 contain the values of channel 1 and channel 2 after they 
have been sorted out in DAT. They should be dimensioned to dim(DAT)/NCHAN. 

0FL01 and 0FL02 are loaded at the end of a record and the beginning 
of the next record if the end of the record falls within the integration 
domain. In order to conduct all the operations required in these overflow 
arrays, 0FL01 must be dimensioned to at least 2 X INTREQ + IBIREQ + 1, and 
0FL02 must be at least as large as INTREQ + IAP2PI/2 + 2 X IAGREQ + IBIREQ + 2. 

2. Integration' Sequences 

In order to process the data as quickly as possible, an extremely 
low order integration process has been used. In fact the method used is 
[actually simply the sunning of a histogram. It is felt that this does not 
add appreciably to the error of the final output, however, due to the re- 
latively small number of sample points in any image which the frequency 
response of the original data recording process allows. If a much slower 
scan rate or a much higher frequency response were possible, then a higher 
order integration process, for example any of the Newton Cotes methods 
[Ref. 18], might be used. 

The integration process linearizes the reference frame in which 
beam wander and image shape are represented by including the factor 
C2(I+IDX), where IDX is approximately equal to a quarter period of C2. That 
is, letting x represent C2 and y represent Cl, 
x = sincot = simrl/IPI. 

dx/di = (tt/IPD-COSttI/IPI = (ir/IPI ) * si n (tt 1/ IPI+IP 1/2 ) 

= (ir/IPI)-x(I+IDX). 

Thus, 

ZC1(I)C2(I)C2(I+IDX) - /x-y-(dx/dl)-dl . /x-y dx 
ZC1 (I )C2( I+IDX) " /y-(dx/dl)-dl fy dx 
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3. CPU Time and Output 

WANDER/SPREAD requires about 1 minute 30 seconds of CPU time and 
outputs approximately 350 punched cards for a 15 second interval of data 
collection with a scan rate of 120 images/second. 
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D. GRAPHICAL OUTPUT 

Given below is BEAD, a program which plots time average image shape with 
and without beam wander, and three intervals of beam wander. The duration 
of these intervals is one second with the program configured as shown. The 
variables WCAL and SCAL scale the values of beam wander and the abscissa of 
the image plot. They are determined from the fraction of the scan drive 
signal in the distance between diffraction maxima generated by the cali- 
bration grid described above, and from the amplitude of channel 2 taken from 
either the printed output of WANDER/SPREAD or 7T09. TCAL is a factor to 
scale the graph of- beam wander to 1 second. 

This program is based almost entirely on Ref. 19, which explains the 
various subroutines called by the program. Figures 13, 14, and 15 show 
typical graphs plotted by this program. 

E. FREQUENCY ANALYSIS 

The program, BEAT, shown below, takes the Fourier transform of the beam 
wander, using the library subroutine RHARM. RHARM is not precompiled and 
must be requested in the Job Control Statements submitted with the program. 
It is based on the Cooley-Tukey algorithm [Ref. 20], and is suitable for 
processing a set with a large number of data points. Figure 16 shows a 
typical frequency spectrum plotted by the program. 

F. INSTRUMENTAL AND COMPUTATIONAL ERROR 

As can be seen from Figure 13, a strong 60 cycle jitter is introduced 
by either the apparatus or the method of data processing. Which of these 
factors is the primary cause is not known, since all the data was taken 
using the scan rate of 60 Hz. In either case the magnitude of this jitter 
is a convenient measure of the overall data collection and reduction process 
error. 
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FIGURE 13. Typical Plot of One-Second of Beam Wander 




FIGURE 14. Typical Image Shape. Half-Width = 2 Microradians 
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FIGURE 15. Typical Image Shape Without Beam Wander Removed. 

Half-Width = 2 Microradians 
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FIGURE 16. Typical Frequency Spectrum of Beam Wander 



In addition, a strong component of the beam wander spectrum at about 
20 Hz can be seen in Figure 16. This frequency corresponds to the length 
of the record and sample rate used in the digitizing process. It is be- 
lieved to be entirely the result of data reduction, since it appeared in 
practically every file processed, always corresponding to the sampling rate 
and record length in use. This could probably be eliminated in the WANDER/ 
SPREAD program by recomputing the bias of C2 at every image. However, this 
would likely increase the magnitude of the 60 cycle jitter and would be 
extremely costly in CPU time. 

Figure 17 shows a plot of beam wander taken using alternate images. As 
can be seen the 60 cycle jitter is substantially reduced. This indicates 
that the jitter is at least partly the result of non-uniformity in the scan 
between the forward and reverse directions. 

In general it has been found that better results have been achieved in 
the data processing if the scan amplitude was small. This allowed more 
sample points within each image for a given digitizing rate, and thus a 
better statistical base for data reduction. This is particularly true in 
the case of computing image shapes. 

Table I shows the RMS beam wander values for the data that was processed. 
The image shapes were not gaussian in most cases, and so are not included. 
Where the shapes roughly approximated a gaussian distribution, enlarged 
versions of Figure 14 were used to obtain the half width of the image at 
its half maximum point. These varied from 2.0 to 11.2 microradians with an 
average of 4.0 microradians for the eight files for which this was possible. 
The image shapes with beam wander included varied from 1.5 to greater than 
55.0 microradians for six files which were roughly gaussian. 

The fact that many of the image shapes appeared skewed is believed 
primarily due to the receiver not having been located in the center of the 
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FIGURE 17. One Second of Beam Wander 
Using Alternate Images 



gaussian distribution. Hallway trials demonstrated the necessity for the 
telescope being centered. However, this was not always feasible in the 
field due to the propagation distance and the fact that all the experiments 
were conducted in conjunction with the scintillation work reported in 
Ref. 13, the apparatus for which also needed to be near the center of the 
beam. 



73 



TABLE I 



Date 


Time 


[Ref. 13] 


RMS beam 
Wander 


Image half 
B.W. Removed 


width 
Not Remc 


9 Oct 


2049 


- 


4.6 


2.0 


2.0 




2153 


0.307 


6.1 


2.9 


1.5 




2216 


0.294 


7.8 


- 


- 




2255 


0.271 


7.8 


7.7 


- 


18 Oct 


2020 


0.292 


18.5 


- 


- 




2022 


0.292 


20.0 


- 


- 




2044 


0.375 


14.5 


4.4 


- 




2049 


0.407 


13.5 


4.3 


- 




2122 


0.503 


19.5 


5.3 


18.2 




2151 


- 


24.4 


2.2 


13.8 


22 Oct 


1626 


0.506 


13.0 


11.2 


12.8 




1636 


0.533 


12.2 


- 


- 




1649 


0.555 


30.2 


- 


- 




1710 


0.510 


13.9 


- 


>55 




1718 


0.533 


11.1 


- 


- 




1729 


0.563 


27.4 


- 


- 


All beam wander and 


half width 


values in microdians. 




All C 


values in 10* 


■7 -1/3 

m 
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V. CONCLUSIONS 



The system described above successfully operated to its design criteria 
during the three trials. Consistent measurements of beam wander values in 
the range of 4.6 to 30.2 microradians and image half widths from 2.0 to 11.2 
microradians testify strongly for the resolution of the apparatus. 

Thus, there seems no reason why continued development of the system 
should not go forward with the basic concept, i.e., measurement with a 
scanning telescope, as a basis. 

The experiences related above indicate that the following areas for such 
development would be the most productive at present. 

The telescope should be reconfigured for simultaneous measurement in the 
visible and infra-red regions. 

An automatic tracking system for the telescope needs to be constructed, 
in order that it may be operated with the shipboard laser stabilization 
system, which is now nearing completion. This is necessary in order to 
obtain measurements which are as free as possible of influences due to 
anomalies in the turbulence conditions brought about by shoreline terrain. 

In this regard, it should be mentioned that conversations with sailing enthus- 
iasts have indicated that conditions over the Monterey Bay propagation path 
are not normally homogeneous, but rather that in about the western two- 
thirds of the path a strong northwesterly breeze usually blows, while in the 
eastern third of the path, the air is nearly dead calm. In addition, to 
provide good correlation of data with C n measurements, a propagation path as 
nearly perpendicular to the prevailing wind as possible is desirable. The 
shipboard transmitter will provide the ability to seek out such conditions. 
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Further improvements in the data processing system would greatly in- 
:rease the volume of information that could practicably be obtained from 
experiments. It has been found that with the system described above only 
about a minute of experimental data can be processed per day, primarily 
because of competition with other users of the IBM 360. An automated 
system using a small computer such as those employed in the experiments 
described in Refs. 14 and 8 would be an excellent solution to this problem. 

Finally, it should be noted that the present status of the project is 
approximately six months ahead of what was planned at the outset. Thus, 
it would appear that, with continued effort, there is no reason why the 
remaining work just described cannot be successfully accomplished. 
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